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Abstract 

Let S\, S2, S3, S4 be given finite sets of pairs of n-by-n complex 
matrices. We describe an algorithm to determine, with finitely many 
computations, whether there is a single unitary matrix U such that 
each pair of matrices in Si is unitarily similar via U, each pair of 
matrices in S2 is unitarily congruent via U, each pair of matrices in 
53 is unitarily similar via U, and each pair of matrices in 54 is unitarily 
congruent via U. 
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1 Introduction 

Our goal is to solve the following 

General Problem. Let <Si, S2, S3, 1S4 be given finite sets of pairs of n-by-n 
complex matrices. Describe an algorithm to determine, with finitely many 
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computations, whether there is a single unitary matrix U such that each pair 
of matrices in Si is unitarily similar via U, each pair of matrices in 5 2 is 
unitarily congruent via U, each pair of matrices in S3 is unitarily similar via 
U, and each pair of matrices in S4 is unitarily congruent via U. 

This General Problem includes as special cases the problem of deter- 
mining whether finitely many pairs of matrices are simultaneously unitar- 
ily similar (<5>2 = S3 = S4 = 0) a.s well as the problem of determining 
whether finitely many pairs of matrices are simultaneously unitarily congru- 
ent (Si = S3 = S4 = 0). 

All of the matrices that we consider are complex and square. Two ma- 
trices A and B of the same size are unitarily similar if there is a unitary 
matrix U such that A = UBU*; they are unitarily congruent if there is 
a unitary matrix U such that A = UBU T . Given pairs of n-by-n matri- 
ces (Ai, Bi), . . . , (A m , B m ) are simultaneously unitarily similar if there is a 
unitary matrix U such that Aj = UBjU* for each j = 1, . . . ,m; they are 
simultaneously unitarily congruent if there is a unitary matrix U such that 
Aj = UBjU T for each j = 1, . . . , m. The trace of a matrix A is denoted by 
tr A We adopt the notation and terminology of [I]. 

2 Unitary similarity of a pair of matrices 

Any finite formal product of nonnegative powers of two noncommuting vari- 
ables s, t 

W(s,t) = s m n ni s m2 t n2 ■■■s m H nk , m 1 ,n 1 ,...,m k ,n k >0 

is a word in s and t. The sum mi + ni + m,2 + n2 + • ■ • + ra k + n k is the length 
of the word W(s,t), and the nonnegative integers rrii and n« are its factor 
exponents. A word in A and A* is 

W(A, A*) = A mi {A*) ni A m2 (A*) n2 ■ ■ ■ A mk {A*) n " (1) 

If A = UBU* for some unitary U £ M n , a calculation reveals that 
W(A, A*) = UW(B,B*)U*, so W(A, A*) is unitarily similar to W(B,B*). 
Thus, unitary similarity of A and B implies that 

trW(A,A*) =trW(B,B*) (2) 

for every word W(s,t) in two noncommuting variables. 
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A theorem of W. Specht [7] provides a converse for this implication. Vari- 
ous authors have provided bounds to show that only finitely many words need 
to be considered [6]; the bound in the following theorem is due to Pappacena. 


Theorem 1 Let complex n-by-n matrices A and B be given. The following 
are equivalent: 

(a) A and B are unitarily similar; 

(b) tr W(A, A*) = tr W(B, B*) for every word W(s, t) in two noncommuting 
variables; 

(c) tr W(A, A*) = tr W(B,B*) for every word W(s,t) in two noncommuting 
variables whose length is at most 

/ 2n 2 1 n 

n \ 7 + 7+ o - 2 - 

V n - 1 4 2 



3 A basic lemma 



The key to obtaining our criteria for simultaneous unitary similarity and si- 
multaneous unitary congruence is understanding the consequences of certain 
intertwining relations involving a special block matrix. 



Lemma 2 Consider the complex k-by-k block matrices 
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in which every block is n-by-n. Define A iii+ i = B iyi+ i = I n for all i = 
1, . . . , k — 1 and A^ = By = whenever % > j. Let W = [W^-]^ - =1 be an 
nk-by-nk matrix that is partitioned conformally to A and B. 

(a) Suppose that AW = WB. Then W is block upper triangular and Wu = 
W 22 = --- = W kk . 

(b) Suppose that W is unitary and AW = WB, that is, A and B are unitarily 
similar and A = WBW*. Then W u = U is unitary and W = U®U®- ■ -@U 
is block diagonal. Moreover, Aij = UBijU* for all i and j. 

(c) Suppose that AW = WB. Then W is block upper triangular, Wu = Wu 
if i is odd, and Wu = Wu if i is even. 

(d) Suppose that W is unitary and AW = WB, that is, A is unitarily con- 
gruent to B and A = WBW T . Then Wu = U is unitary, Wu = U if i is 
odd, Wu — U if i is even, and W = U®U@U®--- is block diagonal. 
Moreover, 

(dl) A^ = U Bij U* if i is odd and j is even; 
(d2) Aij = UBijU T if i and j are both odd; 
(d3) Aij = UBijU* if i and j are both even; and 
(d4) A^ = UBijU T if i is even and j is odd. 

Proof. A computation verifies the assertions in (a) and (c) about block 
triangularity: compare blocks in the respective identities AW = WB and 
AW = WB, starting in block position (k, 1). Work to the right until reaching 
block position (k, k — 1). Move up to block position (k — 1, 1) and work to 
the right until reaching block position (A; — 1, k — 2). Repeat this process, 
moving up one block row at a time, until reaching block position (2, 1). 

The assertions in (a) and (c) about the main diagonal blocks of W follow 
(once one knows that W is block upper triangular) from comparing blocks 
in the respective identities in positions (1, 2),. . . ,(k — 1, A;). 

The assertions in (b) and (d) about W reflect the facts that a block 
triangular unitary matrix is block diagonal, and the direct summands in a 
unitary direct sum are unitary. The asserted relationships between A^ and 
Bij follow (once one knows that W is block diagonal) from the respective 
identities AijWjj = WuAij and AijWjj = WaAij. ■ 

In summary, if the matrices A and B in ((3]) and (jl]) are unitarily similar 
then all of the pairs (Ay, B^) are simultaneously unitarily similar. If A and 
B are unitarily congruent, then there is a single unitary matrix U involved in 
four types of unitary equivalence: certain pairs (A^, B^) are simultaneously 
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unitarily similar via U or U , and certain pairs are simultaneously unitarily 
congruent via U or XJ. 

4 Simultaneous unitary similarity 

It is useful to have an explicit statement of the criterion for simultaneous 
unitary similarity that is implicit in Lemma [2J 

Theorem 3 Let pairs (Ax, Bx), . . . , (A m , B m ) of n-by-n complex matrices be 
given. Choose k large enough so that the matrix A in has at least m blocks 
above the second block superdiagonal. Place the matrices Ax, ... , A m in those 
blocks in any order, and place zero matrices in any unfilled blocks. Place 
Bx, ... , B m and zero matrices in corresponding blocks of the matrix B in (TJJ). 
Then A is unitarily similar to B, if and only if the pairs (Ax, Bx), . . . , (A m , B m ) 
are simultaneously unitarily similar. 

For example, one could choose k = m+2 and place the respective matrices 
of the pairs (Ax, Bx), . . . , (A m , B m ) in positions (1,3),..., (1, k) of the first 
block rows of © and ®, or in the blocks (1,3), (2,4), . . . ,(k — 2,k) of the 
third block superdiagonal of and (jl]). 

A natural extension of Specht's criterion to more than a single pair of 
matrices follows from the preceding theorem. 

Corollary 4 Given pairs (Ax, Bx), . . . , (A m , B m ) of n-by-n complex matrices 
are simultaneously unitarily similar if and only if tr w(Ax, A*, . . . , A m , A^) = 
tvw(Bx, B\*, . . . , B m , B^) for all words w(sx, tx, ■ ■ ■ , s m , t m ) in 2m noncom- 
muting variables. 

Proof. Let k = m + 2 and consider a matrix A of the form ([3]) that is 
constructed by placing the matrices Ax, ... , A m sequentially in the m blocks 
of its third block superdiagonal, and placing zero blocks in all of its other 
blocks above the third block superdiagonal. Construct a matrix B of the 
form (j3J) in the same way using the matrices Bx, ... , B m . 

Consider the following assertions: 

(a) trw(Ax,A* 1 ,...,A m ,A*J = tr w(Bx, B{, . . . , B m , B* m ) for all words 
w(sx, tx, ■ ■ ■ , s m , t m ) in 2m noncommuting variables; 

(b) tr A*) = tr W(B, B*) for every word W(s, t) in two noncommuting 
variables; 
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(c) A is unitarily similar to B; 

(d) The pairs (Ax, Bx), . . . , {A m , B m ) are simultaneously unitarily similar. 
It suffices to show that these four assertions are equivalent. 

(a) =>■ (b) Each block of any word W(A,A*) is a linear combination of the 
identity matrix (we may think of it as an empty word) and words of the form 
w(Ax, A\, — , A m , A m ). 

(b) =>- (c) Theorem [TJ 

(c) =>- (d) Theorem [3j 

(d) =>• (a) The same computation that verified the identities ([2]). ■ 

Simultaneous unitary similarity of given pairs (Ax, Bi), . . . , (A m , B m ) of 
ro-by-n complex matrices is equivalent to unitary similarity of two particular 
block matrices; Theorem [TJ ensures that this latter unitary similarity can 
be confirmed or refuted with finitely many computations. Thus, there is a 
finite algorithm to determine whether a finite number of pairs of matrices 
are simultaneously unitarily similar. 

Likewise, the trace criterion in Corollary [U requires only finitely many 
computations: it requires verification only of enough identities of the form 
txw(Ax, A\, . . . , Am, A m ) = trw(Bx, B*, ... , B m , B m ) to ensure satisfaction 
of the finite number of identities of the form tr W(A, A*) = tr W(B,B*) 
required by the Pappacena upper bound in Theorem [TJ 

5 Unitary congruence of a pair of matrices 

Before we consider the role of Lemma [2] in assessing simultaneous unitary 
congruence and other simultaneous unitary equivalences, we need to consider 
the simplest case of unitary congruence of a single pair of matrices. If U is 
unitary and A = UBU T , a calculation reveals that 

AA* = U(BB*)U*, AA = U(BB)U*, and A T A = U(B T B)U*, (5) 

so three pairs of matrices related to A and B are simultaneously unitarily 
similar. Fortunately, this necessary condition is also sufficient: 

Theorem 5 Complex n-by-n matrices A and B are unitarily congruent if 
and only if the three pairs (AA* , BB*), (AA,BB), and (A T A,B T B) are 
simultaneously unitarily similar. If either A or B is nonsingular, the third 
pair may be omitted. 
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Proof. See 0. ■ 

The preceding theorem and Theorem [3] imply the following criterion for 
unitary congruence. 

Theorem 6 Complex n-by-n matrices A and B are unitarily congruent if 
and only if the An-by-An matrices 
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are unitarily similar. If A or B is nonsingular, they are unitarily congruent 
if and only if 



K' 



are unitarily similar. 



In 


AA* 


AA 





In 










In 










and K'n 



/„ BB* BB 









In 





(7) 



When one applies the criterion in Theorem [HI the bound in Theorem [T] 
ensures that it suffices to verify identities of the form 



trW(K A ,K* A ) = tr W(K B ,K* B ) 
for all words W(s,t) of length at most 



An 



32n 2 
An - 



- + 2n-2. 
A 



(8) 



(9) 



The matrices ([6]) and (J7]) are nilpotent of index four, so in verifying the trace 
identities (IE|) we need to consider only words, all of whose factor exponents 
are three or less. 

For n = 2, the bound says that it suffices to consider all words of 
length at most 37, so a great many words must be considered, even in the 
smallest case. For n = 3 the upper bound on the length is 66; for n = A it is 
116. In contrast, when employing Specht's criterion to test a pair of 2-by-2 
matrices for unitary similarity, it suffices to check traces of only three words 
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of length at most two; to test a pair of 3-by-3 matrices it suffices to check 
traces of only seven words of length at most six. It is not known whether 
the special form of the matrices in (j6]) and (J7|), some special features in low 
dimensional cases, or some clever insight would permit the upper bound (jUJ) 
to be reduced significantly. 



6 Simultaneous unitary equivalences 

We now make explicit the solution of the General Problem that is implicit 
in Lemma 121 

Theorem 7 Let m\, m^, m^, and 1714 be given nonnegative integers. Let 

( aW d(1)\ (a(X) d(1)\ ( a(2) r (2)x ( A (2) r (2)x 

K^-l j D \ )•> ■ ■ ■ A 71 "*!! D m\)i \^\ , £>i )i ■ ■ ■ i \^mi-, £>m2ji 

(Ai , B x ^ ) , . . . , (A m \ > B m \ ) , (A 1 \ b[ ), . . . , (A m \ , B m \ ) 

be given pairs of n-by-n complex matrices. Choose k large enough so that the 
matrix A in (T3|) has enough blocks above the second block superdiagonal to 
accommodate the following construction: 

(1) Place the matrices A x , . . . , Ami (i n an V desired order) in (i,j) blocks of 
A such that i is odd, j is even, and j — i > 2; place the matrices B^\ . . . , Bml 
in corresponding positions in B. 

(2) Place Af\...,A^l (in any desired order) in (i,j) blocks of A such that 

(2) (2) 

i and j are both odd and j — i > 2; place B\ ,.. . , Bmi in corresponding 
positions in B. 

(3) Place Af\.. .,A$ 3 (in any desired order) in (i,j) blocks of A such that 

(3) (3) 

i and j are both even and j — i > 2; place B\ , . . . , B m ^ in corresponding 
positions in B. 

(4) Place A?\...,a!& (in any desired order) in (i,j) blocks of A such thati 
is even, j is odd, and j — i > 2; place b[ 4 \ . . . , B^ in corresponding positions 
in B. 

(5) Place zero matrices in any unfilled blocks of A and B. 

Then A is unitarily congruent to B if and only if there is a unitary matrix 
U such that 

(V) Af ] = UB\ l) U* for alli = l,...,m x ; 
(2') Af ] = UBf ] U T for alli = l,...,m 2 ; 
(3') Af ] = UBf ] U* for alii = 1,..., m 3 ; and 
(A') Af ] = UB^ ] U T for alii = l,...,m 4 . 
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Suppose that four sets of pairs of n-by-n complex matrices are given 
(some of these sets may be empty), and it is required to determine if the 
simultaneous unitary equivalences stated in the General Problem are valid. 
Our algorithm proceeds as follows: 



Algorithm 8 Construct the block matrices A and B according to the pre- 
scription in the preceding theorem. Then construct 



K, 



I n AA* 
/„ 



AA 
A T A 

In 





and K B 



I„ BB* 



L 



BB 
B T B 

In 





The four given sets of pairs of matrices satisfy the required simultaneous 
unitary equivalences if and only if Ka and Kb are unitarily similar. That 
unitary similarity can be confirmed or refuted with finitely many computations 
by using Theorem^ In those computations, only words with factor exponents 
at most 3 need to be considered. 



7 Some comments on previous work 

A criterion for simultaneous unitary similarity that reduces the problem to 
one of verifying unitary similarity of a single pair of block matrices is in Sec- 
tion 2.3 of the 1998 paper [Sj. We employ different block matrices in Lemma 
[2] because we want it to embrace both simultaneous unitary congruence and 
simultaneous unitary similarity. 

The problem of simultaneous unitary similarity was also studied in the 
2003 paper [T]. The criterion developed there is formally equivalent to the 
finite version of our Corollary HI 

The recent paper [2] makes the important observation that the criterion 
in [3] can be combined with a test for simultaneous unitary similarity to give 
a criterion for unitary congruence that can be verified with finitely many 
computations; the authors use the test in [TJ. However, because of an extra 
condition imposed in [1] on matrix families whose simultaneous unitary simi- 
larity is to be tested, the criterion in [2] for determining unitary congruence of 
two general matrices requires assessing simultaneous unitary similarity of the 
four pairs (AA*, BB*), (AA, BB), (A T A* , B T B*), and (A T A, B T B) rather 
than the three pairs in (jSj). 
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